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Abstract 

The Proca wave equation describes a classical massive spin 1 particle. We 

analyze the gravitational interaction of this vector field. In particular, the 

spherically symmetric solutions of the Einstein-Proca coupled system are 

obtained numerically. Although at infinity the metric field approaches the 

usual Schwarzschild (Reissner-Nordstrom) limit, we demonstrate the absence 

of black hole type configurations. 
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I. INTRODUCTION 



A massive vector meson (spin 1 particle with a non-trivial mass) is described by a one- 
form field which obeys the Proca wave equation ffl. Early development of the Proca the- 
ory was concerned with the classical and quantum electrodynamics of a massive photon. 
However, strong experimental limits on photon's mass (see, e.g., 0) in combination with 
theoretical arguments based on the idea of gauge invariance (which ultimately led to the 
standard model of electroweak interactions) have closed the electrodynamical chapter in 
the history of this theory. A further discussion of the differences between the Proca and 
electromagnetic fields can be found in Q. 

At present, interest in the Proca field is twofold. Firstly, the Proca model presents a 
convenient theoretical "laboratory" for the study of Lagrangian and Hamiltonian theories 
with second class constraints 0. Secondly, although it is irrelevant for the electrodynamics, 
a massive vector meson often appears in the spectra of many non-trivial field theoretical 
models, including some classes of generalized theories of gravity. In connection with this, 
it is interesting to investigate the specific physical effects arising in such models due to the 
interaction of Proca particles with electromagnetic, gravitational and other physical and 
geometrical fields. 

The interaction of spin 1 field with electromagnetic field is known to be free of algebraic 
inconsistencies as well as of acausal wave propagation (v > c) when the coupling is minimal 
(or modified by the addition of an anomalous magnetic dipole moment). However, acausal 
propagation anomalies arise for more general interaction Lagrangians || . Similarly, acausal 
propagation takes place (along with algebraic inconsistencies) for a Proca field coupled 
minimally to external torsion field ||. Different aspects of the interaction of classical and 
quantum vector field with torsion have been analyzed recently in J7|. 

Early studies of the gravitational interaction of the Proca field were centered around the 
black hole issue. Qualitative analysis of the self-consistent Einstein-Proca system revealed 
the absence of black hole type solutions (possessing a regular horizon) with an external vec- 
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tor meson "hair" PHlO|j. At the same time several exact spherically symmetric solutions of 
the Proca wave equation on the classical Schwarzschild background spacetime were obtained 



I 1 ^L3[ . Assuming that a massive vector field source was located on a thin spherical shell 
outside the Schwarzschild horizon, it was demonstrated in [|11[] that the meson field may 
change the structure of the spacetime near the central singularity. For point vector field 
sources located at the origin [fT^l , or at a finite distance from the origin [13|, it was shown 



that the range of the meson field is reduced by the metric gravitational field. The energy- 
momentum invariant was found to be divergent on the Schwarzschild horizon. However, it 
should be noted that, contrary to the Abelian Proca case, the non-Abelian massive vector 
field (with mass of a Yang-Mills field coming from a spontaneous symmetry breaking mech- 



anism) may form a black hole type configuration |1J|. The results of numerical analysis of 
the spherically symmetric gravitationally interacting complex spin 1 field have been reported 
recently in |]I5 |. In this case the Einstein-Proca system admits everywhere regular "boson 



star" -type solutions (cf. with massive scalar boson stars |T6|jri 

A direct motivation for our current study comes from the metric-affine theory of gravity 
(MAG). In Einstein's general relativity the spacetime geometry is described by the curva- 
ture 2- form R a " . In MAG two post-Riemannian structures are introduced: the 1-form of 
nonmetricity Q a p and the torsion 2-form T a . For a comprehensive review of this theory 
sec 



18 1 . Already the early investigations [0,01 °f the models with the simplest possible 



MAG Lagrangians, which include only a linear Hilbert term, quadratic segmental curvature 
invariant, and a single trace torsion or Weyl nonmetricity square term, have shown that an ef- 
fective Einstein-Proca theory arises naturally from the vacuum MAG field equations (cf. also 
fl2| , |2T|| ) . This result was subsequently extended to a very general family of MAG Lagrangians 
|22|-f2~i||. In all these models the effective Proca field describes the triplet of post-Riemannian 
one-forms which are proportional to each other: the Weyl covector Q := g^Qap/^, the tor- 
sion trace T := e a \T a , and the nonmetricity one- form A := ^ a e^\Q a/ 3 — Q. The mass of the 
effective vector particle is constructed from the coupling constants of the MAG Lagrangian. 



For a complete review of the known exact solutions of MAG see |25 



In this paper we study the spherically symmetric static solutions of the coupled Einstein- 
Proca system of field equations. A preliminary analysis of the limiting cases of this problem 
shows a possibility of solutions which combine the exponential "Yukawa" type behavior of 
the Proca potential at the origin with the asymptotically Schwarzschild solution far away 
from the source. We will present the corresponding solutions which have been obtained by 
the application of numerical integration techniques. 



Our main conventions and notation are taken from [fig] . In particular, the 77-basis of 
the exterior algebra is constructed from a coframe one-form $ a with the help of the Hodge 
duality operator: <q a = ^ = *($ a A ^),7] a ^< = A# A $ 7 ). The dual frame is 
denoted as e a . The Greek indices a, ft, . . . = 0, . . . , 3 label anholonomic components, and 
the metric signature is (—,+,+, +). 



II. EINSTEIN-PROCA THEORY 

The Lagrangian four-form of the Einstein-Proca system reads 

V = - — R a(3 A Va s - - (dA A*dA + m 2 AA*A) , (2.1) 
2k 2 ^ > 

where k is the gravitational constant (k = £ 2 ) and m is the rest mass of the vector field A. 
The corresponding field equations arise from the independent variation of the action with 
respect to the coframe and the Proca one-forms, and read: 

d*dA + m 2 *A = 0, (2.2) 
^R* AT tafh = KE a . (2.3) 

Here the canonical energy-momentum three-form of the massive vector field 

S a = - {(e a \dA) A *dA - {e a \*dA) AdA + m 2 [(e a \ A)*A + (e a \*A) A A]} , (2.4) 

represents the usual source of the gravitational field. 

It is worthwhile to recall the relationship of ( |2.1| )-( ]2TB| ) to the MAG theory. As we 
have mentioned already, the same physical system arises in MAG as an effective system 
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T9| , |20|j22| -p5| in which the effective covector Proca field is (in the notations of our previous 



paper 

A = ^k <j). (2.5) 

Here <ft determines the three nontrivial post-Riemannian pieces of nonmetricity and torsion 
(the triplet of one-forms) 

(l) T a = (3) T a = 0> (l)g^ = (2) Qa/3 = ^ ^ 

Q = k (f>, A = h(f>, T = k 2 <p. (2.7) 

The effective mass m 2 of the vector particle and the constants ko,ki,k 2 are constructed 
from the original coupling constants of the MAG Lagrangian which contains all possible 
quadratic invariants of the torsion and nonmetricity (11 terms) together with the linear 
Hilbert type term (multiplied by the constant k) and the Weyl segmental curvature quadratic 
term (multiplied by the constant z±). See []23| and for more details (note, however, that 



in the present paper we assume that the cosmological constant is zero) 



III. SPHERICALLY SYMMETRIC STATIC CASE 

In terms of the local time and space coordinates (r, r, 9, 0), the general spherically sym- 
metric ansatz for the coframe can be written as 

$o = fd T , d l =jdr, § 2 =rd6 ) = r sin0 d(j) . (3.1) 

The geometrical meaning of the function g(r) becomes evident when one computes the 
volume four-form 

rj = $ 6 A ^ A d 2 A 9 = g r 2 sin 6dr A dr A d6 A d<p. (3.2) 

Thus, g(r) measures the deviation of rj from the standard spherically symmetric spacetime 
volume form. In a regular oriented spacetime domain we naturally have to assume 
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< g(r) < oo. 



(3.3) 



The general static spherically symmetric configuration of the coupled Einstein-Proca 
system is described by the three functions / = f(r),g = g(r), and u = u(r) which enter the 
spherically symmetric ansatz for the Proca field as follows 



A = — $° = - dr. 

rf r 



Substitution of (|3.1|)-(|3T4|) into the Proca field equation fl2.2|) results in 



11 

r 2 g 



2 , v / 



9 \ r 



' 2 ' 

m u 



rf 



or, equivalent ly 



// 9' ( , u\ m 2 g 2 
u I u ) ^— u — 0. 

r 



g \ rj f 2 
A direct calculation of the energy-momentum 3-form yields 



1 



2r 2 g 2 f 2 
1 



-f[u'--) +m 2 g 2 u 2 \Sl m 



+ 



2rV/ 2 

For the sake of completeness we also write down the Einstein 3-form 



~ R* A rj a/3j = 

2 r^g 1 



f g 2 ' 

2 ~ r + l ~T 2 



f 



P 



r«g« J + ( r / ) ~ ^ J-g + 2r 7 " r * 1 @ % + ^ ' (3 ' 8) 



-2„2 



r 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



Inserting ( |3.7| ) and (|3.8|) into (|2.3|) we find (after some algebra) that the equations cor- 
responding to a = 6, 1 are 



2 (r (/ 2 )' + / 2 - g 2 ) = k 



I 2 

f 2 ' 



29 2 

m —u — k\u 



u 



{ 2\l 29 2 

r (g ) = nm —u . 



(3.9) 
(3.10) 



Furthermore, it is straightforward to see that the (second order) equation corresponding to 
a = 2, 3 is a consequence of ( |3 . 9|) - ( pTTQ|) and 



The scalar invariant |E| := A characterizes the "magnitude" of the energy- 

momentum of the massive vector field. Using (|3.7|), we find that 



IV. PRELIMINARY ANALYSIS 

Before we start the study of the complete system it is instructive to recall two particular 
cases: namely, massless vector particle in curved spacetime and massive vector particle in 
Minkowski spacetime. 

For the massless vector particle 

m = 0, (4.1) 

and one finds, from (|3.10| ), that g — go — const. Consequently, the vector field equation can 
now be easily integrated to give the usual Coulomb solution 

u = q = const => A = - dr. (4.2) 

r 

[Strictly speaking, the general solution reads u = q + (3r, but one can put (3 = since it 
contributes only an exact form to A]. Turning to (|3.9| ), we immediately recover the well 
known Reissner-Nordstrom solution 

M 1 -"+^)* (4 - 3) 

with the integration constant M interpreted as the mass of the gravitating source and q/go 
as its electric charge. 

On the other hand, for the Minkowski spacetime the metric functions are / = g = 1 and 
the solution of the Proca field equation (|3.6|) yields the well known Yukawa potential 



-mr 

u = q e - mr A=- dr. (4.4) 
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This shows that the vector field is practically zero at distances much greater than the typical 
length r = 1/m. 

We expect that a spherically symmetric configuration of coupled Einstein and Proca 
fields will combine both of the typical features of the above limiting cases. Namely, for small 
values of mass m there will be a large part of space inside the sphere of radius r = 1/m 
where the function u is to a high degree of approximation constant. In this region the exact 
solution will naturally be approximated by ft4.2|) and the metric will assume the familiar 
Reissner-Nordstrom form ft4.3| ). However, due to its massiveness the field A will remain, also 
in curved spacetime, confined to a finite spatial volume, whereas for r — » oo one expects a 
fast decay u — > which leaves one with pure Schwarzschild metric. 

The following observation will be very useful in the discussion of exact solutions. Multi- 
plying ( |3.6|) by u/g and using the Leibniz rule one finds that 

where the right-hand side is positive definite in a regular spacetime region. 
Identity (|4.5| ) represents a particular case of the general relation 

db = —{dA A *dA + m 2 A A *A), (4.6) 

where b := —A A *dA. The latter identity holds true for all solutions of the Proca field 
equation ( |2.2| ). 



V. DIMENSIONLESS SYSTEMS 

The general system is nonlinear and apparently cannot be integrated analytically. Con- 
sequently, we will present the results of numerical integration in the remainder of this paper. 
Before we begin the numerical analysis, we introduce a new dimensionless radial variable 

p:=^=, (5.1) 



which allows us to rewrite the system (|3.9|) , (|3.10|) , and ( |3.6|) in the form 
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i>% + *- a )-*T*-&-$- (5 - 2) 

P f p =K%A (5.3) 
d 2 u T ^G 1 dG ( du u\ . A . 

-d? =K F u+ 2Gir P {T p -- p )- (5 - 4) 

Here we have introduced a dimensionless constant 

K := nm 2 , (5.5) 

and defined the functions 

F:=f, G:=g 2 . (5.6) 

Evidently, the dimensionless radial coordinate measures distance from the origin in units of 
the Planck length (k = £ 2 ). At the same time, the parameter \^K, being the ratio of the 
Planck length to the Compton length of the vector particle, characterizes the size of the 
domain where the influence of the Proca field on the spacetime geometry is significant. 

One can, alternatively, study a different dimensionless system after defining the scaled 
metric functions 

F := F/K, G := G/K, (5.7) 

and introducing a new radial coordinate 

£ := y/~Kp = mr = — . (5-8) 
r 

The system (|5.2| )-( [5T4|) then reads 

^H-Wari)'- (59) 

£ — = ^u 2 , 5.10 
d£ f 2 

d 2 u G 1 dG ( du u\ 

ie = 7 u + wii\TCV- (5 ' n) 
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In this form the equations no longer contain a free parameter (such as K) and the new 
dimensionless coordinate measures distance in units of the characteristic ("Compton wave- 
length") scale tq. It is convenient to use both dimensionless systems. The advantage of 
( |5.9p - (|5.1ip lies in the absence of K, whereas the equations ( |5.2| )-( |5^ ) are more transparent 



from the physical point of view when one considers limits of small and big mass m. 

VI. CONDITIONS AT THE ORIGIN AND AT INFINITY 

Before one can start the numerical integration, an appropriate set of initial conditions 
must be specified. Unfortunately, the solution cannot be represented by analytic power 
series expansion for u, F, G at the origin in view of the apparent singularity at p = 0. 

Instead, one can verify that for small values of p, irrespective of the value of K, there is 
an approximate solution of the form 



u ~ q + bp, (6.1) 
2 p 2 

G « yA-J. ( 6 - 3 ) 



F~ (6.2) 



where q, b, c are parameters which determine the initial conditions in the neighborhood of 
the origin p = 0. 



At infinity, p — > oo, following the physical discussion in Sect. [TV], we expect approximate 
behavior of the form 

u — > u exp(— \J~Kp), (6.4) 
F — 9 2 f 1 - — J , (6.5) 
G — gl (6.6) 

where «o,9o are constants. The condition (|6.4|) means that the massive vector field is non- 
trivial only inside a sphere of a finite radius 1/VK ( "Yukawa- type" behavior). On the 
other hand, the conditions ( |6.5|) - (|6.6|) specify purely Schwarzschild asymptotic metric. A 
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more precise form of the limit ( |6.6| ) is easily obtained after substituting ( |6.4|) and ( |6.5|) into 

G « g 2 (l + Ku 2 Ei(-2VKp)) , (6.7) 

x t 

where Ei(rr) = / ^rdt is the integral exponential function. It is worthwhile to recall that 

— oo 

asymptotically, for x — > oo, one has Ei(— x) ~ — e~ x fx. 

We will use the asymptotic conditions (|6.4j) - (|6.6|) , (|6.7|) in the numerical analysis of the 
problem under consideration. 



VII. ABSENCE OF SOLUTIONS WITH HORIZONS 

In this section we show that the spherically symmetric Einstein-Proca system does not 
admit asymptotically flat solutions with horizons. The absence of black holes for a massive 
vector field was first demonstrated by Bekenstein ||. 

Let us consider an arbitrary regular solution uir) which vanishes at two points 7*1 and 
r 2 > r\. u{r\) = u(r 2 ) = 0. Then u(r) = for all r\ < r < r 2 . Indeed, integrating the 
identity ( |4.5[ ) from r\ to r 2 , one finds that 

7 1 ( ( U ' - ») 2 + „» 4 1 * = 14 f«'M - ^] - 14 {An) - ^) = 0. 

(7.1) 

Since the integrand is positive definite, the vanishing of the integral leads to the above 
conclusion. 

Consequently, a nontrivial solution u(r) which vanishes asymptotically at r 2 = 00 (thus 



satisfying the condition ( |6.4| )) cannot have zeros at any finite r\ (since then the solution 
would be trivial: u(r) = for r > r{). 

This leads to the absence of the black hole type solutions of the system (p.2|)-(pT4|). 
In order to see this, let us recall that a black hole necessarily possesses a horizon. Quite 
generally, on a spacetime manifold Ai a horizon is defined as a hypersurface S := {x l G 
M. I a{x l ) = 0} such that: (i) the normal vector rij := <9jO" is null 
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^71^ = 0, (7.2) 

and (ii) S is not an essential singularity. The latter means that all the curvature invariants 
as well as the volume 4-form r\ are nonsingular on the horizon. In particular, the regularity 
of rj follows from the condition ( |3.3| ) on the function g. 

For a spherically symmetric gravitational field configuration, horizon S is evidently a 
sphere a = r = r^. Normal vector is then rii = 8}. Substituting (|3.1| ) into (|7.2j), one obtains 

f 2 (r h ) _ F{r h ) _ 

g\r h ) ~ G(r h ) ~ °" (7 ' 3) 
Since G(rh) is finite in view of fl3.3|) , we find that F must vanish on the horizon S 

F(r h ) = 0. (7.4) 

The last equation formally defines the position of a horizon in the general spherically sym- 
metric spacetime fl3.1|) . Now recall the second requirement: a hypersurface S must be free 
of physical singularities in order to be a horizon. Clearly, the energy-momentum invariant 
scalar fl3.11| ) is regular at r = if and only if 

u{r h ) = 0. (7.5) 

Furthermore, if ( |7.5|) did not hold then the quadratic curvature invariants (obtained by using 
the Einstein field equations (|2.3|) in the definition of |E|) would diverge at because of the 
last term in ( 3.11| ) and (7.4). 



Now we are in a position to conclude that there are no solutions with a horizon and a 
nontrivial massive vector field. Indeed, assume the contrary is true. Then outside a horizon 
S the function u is necessarily given by u{r) = 0, < r < oo, because u vanishes at infinity 
( p.4| ) and at the horizon ( |7.5| ). Consequently, outside S, the system ( |5.9| )-( |5T0| ) has the 



usual Schwarzschild solution G = 1, F = 1 — r^/r. Integrating ( |5.9| )-( 5~TT| ) from to with 
the initial conditions ( |7.4| ) and ( |7.5| ), we find u(r) = everywhere. 

Bekenstein's original proof || was based on the assumption that the three-form b = b a r] a 
defined in (O) is bounded on the horizon. It is easy to see that b a = \ u' — 8° 
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diverges on the horizon ( |7.4| ) unless u vanishes. Thus, for a massive vector field, the form b 
is not only bounded but, in fact, trivial on the horizon. 



VIII. NUMERICAL SOLUTIONS 

After fixing the value of the parameter K to the square of the ratio of the Planck length 
to the Compton wavelength of the vector field, one can start numerical integration at an 
arbitrarily small p with the initial conditions defined by ( |6.1|) -( |6T3"1) . One is free to choose any 
initial value for the "boson charge" function «(0) = q (7^ 0, otherwise u is trivial everywhere). 
Solutions with the correct asymptotic behavior ( |6.4|) - (|6.6|) exist only for fixed values of the 
parameters h = u'(0),c = G(0). Technically, the numerical integration can start at a point 
arbitrarily close to the origin for every chosen values of K and q. In order to obtain the 



asymptotic behavior ( |6.4| )-( p^6| ), a fine tuning of b and c is required which can be achieved 
similarly to the construction of the Bartnik-McKinnon solutions or of the Abrikosov- 
Nielsen-Olesen vortices (see, e.g. and references therein). Alternatively, one can start 
the numerical integration at a sufficiently large radius with the initial conditions taken from 
d6.4D -( fT6l) for arbitrary values of K, M, w . As a cross-check, we have used both integration 
schemes. The resulting approximate solutions turned out to be completely consistent with 
each other. 

Particular solutions for various values of K, q and M are described in Tables H-pH The 
graphical form of the solutions is presented in Figures |],@, and [| In these figures, the 
numerical solutions are depicted for K = 1 and M — 0.1 (dotted lines), M = 0.5 and 
M = 1.5 (broken lines), and M = 2 (solid lines). In all cases we put g = 1 which is always 
possible to achieve by the redefinition of the time coordinate r — > g^r. 

As one can see, the relation between K (formal rest mass of the vector field) and M 
(asymptotic total mass of the solution) plays a decisive role. At the same time, the value of 
the boson charge q at the origin is also important. 



In agreement with the results of Bekenstein et al ]^iTJ| and with the preliminary analysis 
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of Sect. |VII| , all the numerical solutions obtained by us are without horizons. They possess a 
true physical singularity at the origin which provides us with an example of a naked spacetime 
singularity. Stability of these solutions against small perturbations will be studied separately. 

Recalling that the effective Proca field emerges naturally in the general metric-affine 
models, we thus conclude that the presence of the post-Riemannian geometric objects pre- 
vents, in general, a formation of a black hole in MAG theory. Only in the special case when 
the MAG coupling constants are such that the effective mass vanishes, m 2 = 0, the black 



holes can be formed 2C.22 
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FIG. 1. The metric function F = f 2 : K = 1 solutions for the values M = 0.1 (dotted line), 



M = 0.5 and M = 1.5 (broken lines), and M = 2 (solid line). 
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FIG. 2. The metric function G = g 2 : K = 1 solutions for the values M = 0.1 (dotted line), 



M = 0.5 and M = 1.5 (broken lines), and M = 2 (solid line). 
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FIG. 3. The vector field function u: K = 1 solutions for the values M = 0.1 (dotted line), 



M = 0.5 and M = 1.5 (broken lines), and M = 2 (solid line). 
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TABLES 

TABLE I. Solutions with fixed values of q and M 



q = A and M = 0.5 


q = 2 and M = 1 


K 


b 


c 


K 


b 


c 


0.01 


-0.29605 


0.90444 


10~ 6 


-0.00202 


0.99996 


0.10 


-0.66479 


0.69257 


0.01 


-0.16267 


0.84124 


1.00 


-1.20050 


0.34457 


1.00 


-0.56351 


0.08738 


10.00 


-1.67696 


0.08873 


10.00 


-0.69666 


0.01100 



TABLE II. K = 1.00: solutions for q = 1 with different masses M 



Mb c 

0.10 -0.56959 0.76748 

0.50 -0.43469 0.24706 

1.50 -0.13588 0.00528 

2.00 -0.10007 0.00170 



TABLE III. Solutions for M = 1 and different q 





q = 0.50 


q = 1.00 


K 


b 


c 


6 




c 


0.01 


0.60874 


0.08438 


0.25196 




0.34048 


1.00 


-0.01728 


0.00579 


-0.21296 




0.02528 


10.00 


-0.13175 


0.00071 


-0.33056 




0.00308 
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